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We investigate numerically the three-dimensional 0(4) model on 24^-120^ lattices 
as a function of the magnetic field H. We verify explicitly the singularities induced 
by Goldstone modes in the low-temperature phase of the model, and show that they 
are also observed close to the critical temperature. Our results are well described 
by the perturbative form of the model's magnetic equation of state, with coefficients 
determined nonperturbatively from our data. The resulting expression is used to 
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1 Introduction 



The 0{N) spin models (or, more precisely, the 0(A^)-invariant nonlinear cr-models) 
are defined by 

pn = -J s.-s, - H- 5]s,, (1) 

<i,j> i 

where i and j are nearest-neighbour sites on a dimensional hypercubic lattice, 
and Si is an iV-component unit vector at site i. The case N = 1 corresponds to the 
Ising model. We will consider here only > 1. It is convenient to decompose the 
spin vector Sj into longitudinal (parallel to the magnetic field H) and transverse 
components 

Si = S}u + St . (2) 

The order parameter of the system, the magnetization M, is then the expectation 
value of the lattice average S"" of the longitudinal spin components 

M =<1^5j> = <5ll> . (3) 

i 

Two types of susceptibilities are defined. The longitudinal susceptibility is the usual 
derivative of the magnetization, whereas the transverse susceptibility corresponds 
to the fluctuation per component of the lattice average S"*" of the transverse spin 
components 

XL = II = V{< St' > -M') , (4) 



These models are of general interest in condensed matter physics, but have ap- 
plications also in quantum field theory. In particular, the three-dimensional 0(4) 
model is of importance for quantum chromodynamics (QCD) with two degenerate 
light-quark fiavors at finite temperature. If QCD undergoes a second-order chiral 
transition in the continuum limit, it is believed to belong to the same universality 
class as the 3d 0(4) model ^ 0. QCD lattice data have therefore been compared 
to the 0(4) scaling function, determined numerically in ||^. For staggered fermions 
this comparison is at present not conclusive but results for Wilson fermions |^ 
seem to agree quite well with the predictions. 

0{N) models in dimension 2 < d < A are predicted to display singularities 
on the coexistence line T < T^, H = due to the presence of massless Goldstone 
modes In fact, both susceptibilities are predicted to diverge in this region. The 
magnetic equation of state is nevertheless divergence-free, and compatible with these 
singularities. The equation of state was calculated up to order in the e-expansion 
by Brezin et al. 0]. On the basis of this expansion it has been argued that these 
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singularities may be easily observable, since the perturbative coefScient associated 
with the diverging term is quite large. 



A further consequence of the Goldstone singularities is the appearance of strong 
finite-size effects at all T < Tc for H ^ 0. These effects have been studied using 
chiral perturbation theory in ^ . Direct numerical evidence of the Goldstone sin- 



gularities is however lacking, apart from early simulations of the three-dimensional 



0(3) model on small lattices |TT|, where indications of the predicted behaviour were 
found. 

The aim of this paper is to verify explicitly the Goldstone singularities, and to 
investigate their interplay with the critical behaviour and the effect they have on the 
scaling function. We do this by simulating the three-dimensional 0(4) model in the 
presence of an external magnetic field in the low-temperature phase and close to the 
critical temperature Tc. For determining the scaling function, we have also simulated 
at some high-temperature values. First results of our work have been presented at 
Lattice '99 [Q. The plan of the paper is as follows. In the next section we review 
the perturbative predictions for the magnetization and the susceptibilities at low 
temperatures, as well as the analytic results for the magnetic equation of state, 
which is equivalent to the magnetization's scaling function. Our numerical results 
are discussed in Section ^. The fits and parametrization for the scaling function are 
given in Section ^. A summary and our conclusions are presented in Section ^. 



2 Perturbative Predictions and Critical Behaviour 



The continuous symmetry present in the 0{N) spin models gives rise to the so-called 
spin waves: slowly varying (long-wavelength) spin configurations, whose energies 
may be arbitrarily close to the ground-state energy. In two dimensions these modes 
are responsible for the absence of spontaneous magnetization, whereas in > 2 they 
are the massless Goldstone modes associated with the spontaneous breaking of the 
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rotational symmetry for temperatures below the critical value Tc [0. For T < 
the system is in a broken phase, i.e. the magnetization M{T, H) attains a finite value 
M{T, 0) at = 0. To be definite we assume here H > 0. As a consequence the 
transverse susceptibility, which is directly related to the fluctuation of the Goldstone 
modes, diverges as when H —>■ for all T < Tc. This can be seen immediately 
from the identity 

M{T, H) 

XT = . (6) 

This expression is a direct consequence of the 0{N) invariance of the zero-field free 
energy, and can be derived as a Ward identity 0. It is valid for all values of T and 
H. 

A less trivial result pi, M is that also the longitudinal susceptibility is diverging 
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on the coexistence curve for 2 < < 4. The leading term in the perturbative 
expansion for 2 < ci < 4 is if^/^-s 'jj-^g predicted divergence in ci = 3 is thus 



XL 



(T < Tc, H) ~ H-^/^ . (7) 



This is equivalent to an if^/^-behaviour of the magnetization near the coexistence 
curve 

M(T < T„H) = M(T,0) + cH^/^ . (8) 

An interesting question is whether the above expressions still describe the be- 
haviour close to the critical region T < Tc. We recall that the critical behaviour is 
determined by the singular part of the free energy. Its scaling form in the thermo- 
dynamic limit is 

fsit,h) = h-''fs{hy%w-h) . (9) 

Here we have neglected possible dependencies on irrelevant scaling fields and expo- 
nents. The variables t and h are the conveniently normalized reduced temperature 
t = (T — Tc)/Tq and magnetic field h = H/Hq, and 6 is a free length rescaling factor. 
The relevant exponents yt^h specify all the other critical exponents 

yt = ^/i^, yh = ^/i^c; (lo) 

iy, = iy//35, diy = /3il + 6), ^ = /3{5 - 1) . (11) 

Choosing the scale factor b such that b^^h = 1 and using M = —dfs/dH one finds 
the equation 

M = h'/'fcit/h'/^') , (12) 

where fc is a scaling function. It becomes universal after fixing the normalization 
constants Ho and Tq. This scaling function was calculated numerically for the 3d 
0(4) model by Toussaint ||^ and is used in comparison to QCD lattice data 

Alternatively, one may choose = 1. This leads to the Widom- Griffiths form 

y = fix) , (13) 



of the equation of state |14 
where 



y = h/M\ x = t/M^'^. (14) 
It is usual to normalize t and h such that 



/(0) = 1, /(-1) = 0. (15) 

The scaling forms in Eqs. ([12|) and ([T3|) are clearly equivalent. In the following we 
will work with form (|1^) and obtain (|1^) from it parametrically in Section ^. 

The equation of state (|TB|) has been derived by Brezin et al. to order in the 
e-expansion, where e = A—d. Although diverging terms in xt appear at intermediate 
steps of the derivation, they are canceled by diverging xl terms, and the resulting 



3 



expression is divergence-free. This expression has been considered by Wallace and 
Zia in the limit x —>■ —1, i.e. at T < and close to the coexistence curve. In 
this limit the function is inverted to give a; + 1 as a double expansion in powers of 
y and y'^/^-i 

x + 1 = c,y + C2y''/'-' + d,y^ + d^y"/'' + 42/""' + ■ • ■ ■ (16) 

The coefficients Ci, C2 and are then obtained from the general expression of J^. 
The above form is motivated by the if-dependence in the e-expansion of xl at low 
temperatures 0. 

In Fig. |l]we show the function f{x) from 0, its low-temperature [x — *• —1) limit 
and f{x) from [Q, which is obtained from the inverted form of the low-temperature 
expression, Eq. (pIBI) . We see that the low-temperature curve remains close to the 
general one for a significant portion of the phase diagram, including the critical 
point at a; = 0, ?/ = 1 and moving into the high-temperature phase, the region to 
the right of the critical point. The "inverted" curve and the low-temperature curve 
that generated it agree quite well. (We remark that the process of inverting produces 
coefficients that are determined only to order e, while the original low-temperature 
expression was known to order e^.) As mentioned above, the form (^) is equivalent 
to the Goldstone-singularity form for xl at low temperatures if one identifies the 
variable y with the field H. Nevertheless, the fact that this form may describe the 
behaviour also at temperatures close to Tc and higher is not contradictory, since the 
variable y can only be identified with H if M{T, H = 0) 0, which happens only at 
low temperatures. The form (p!6D has explicitly nonnegative values of ?/ at x ~ — 1, 
while the original perturbative expression produces an unphysical negative y in a 
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Figure 1: The function y = f{x) (solid line) from its approximation for x — 1 
(dashed line) and the inverted form (dotted line) from f^. 
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very small neighborhood of this point . Possible problems with the form (|T6[) are 
pointed out in |jl5|, Section 5]. 



As for the large-x limit (corresponding to T > Tc and small H), the expected 



behaviour is given by Griffiths's analyticity condition |jT4 



fix) = Y.^nX^-"^-'^' ■ (17) 



n=l 



None of the curves in Fig. |T] approaches this limit, since the low-temperature curves 
are not valid at large x, and the general curve in its original form is known to have 
problems in this limit 

The perturbative equation of state has been used in p[ to produce pictures of 
the expected behaviour of the 3d 0(4) model for a large range of temperatures and 
magnetic fields. The authors have employed an interpolation of the function from 
1^ with the inverted form from [Q| at low temperatures and the Griffiths condition 
at high temperatures. When compared to Monte Carlo data for the same model in 
1^, the perturbative scaling function shows qualitative agreement. In Section ^ we 
propose a fit of our Monte Carlo data to the perturbative form of the equation of 
state, using (llf 



3 Numerical Results 



Our simulations are done on lattices with linear extensions L = 24, 32, 48, 64, 72, 
96 and 120 using the cluster algorithm of Ref. 0. We compute the magnetization 
M and the susceptibilities Xl,t at fixed J = 1/T (i.e. at fixed T) and varying H. 
We note that, due to the presence of a nonzero field, the magnetization is nonzero 
on finite lattices, contrary to what happens for simulations at H = 0, where one is 
led to consider the approximate form < 1/V^| Sj| >. 

We use the value Jc = 0.93590, obtained in simulations of the zero-field model 
. In Fig. ^ we show our data for the magnetization for low temperatures up to 
Tc plotted versus H^^"^. We have simulated at increasingly larger values of L at fixed 
values of J and H in order to eliminate finite-size effects. The finite-size effects for 
small H do not disappear as one moves away from Tc, but rather increase. 

In Fig. ^ we plot only the results from our largest lattices. The solid lines are fits to 
the form (||), and the filled squares a.t H > denote the last points included in our 
fits. It is evident that the predicted behaviour (linear in H^^"^) holds close to H = 
for all temperatures T < Tc considered. The Goldstone-mode effects are therefore 
observable also rather close to Tc. The straight-line fits coincide with the measured 
points in a wide range of H for low T (for J = 1.2 and 1.1 up to if = 0.32). 
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Figure 2: The magnetization vs. H^^^ in the low-temperature region for fixed J 
1.2, 1.0, 0.95 and Jc and different lattice sizes. 
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Figure 3: The magnetization as a function of H^^^ for fixed J = 1.2, 1.1, 1.0, 0.98, 
0.95 and Jc, starting with the highest curve. The size L is denoted as in Fig. 0. 
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J= 1/T 


M(T, 0) 


Slope c 


if —Range 


xVdof 


0.95 


0.20202(62) 


1.1825(132) 


0.001-0.003 


0.37 


0.98 


0.30650(13) 


0.6495(26) 


0.0005-0.007 


0.38 


1.00 


0.35061(15) 


0.5245(20) 


0.001-0.01 


0.60 


1.10 


0.48257(29) 


0.2940(38) 


0.0015-0.01 


0.54 


1.20 


0.55911(13) 


0.2033(13) 


0.002-0.02 


0.92 



Table 1: Parameters of the fit of the magnetization to M{T, 0) -|- cH^^^. 



With increasing T the coincidence region gets smaller and vanishes at T^. In Table 
m we have listed the fit parameters. The value M(T, 0) obtained from the fits 
is the infinite- volume value of the magnetization on the coexistence line. In the 
neighbourhood of it should show the usual critical behaviour 

M(T < T„ if = 0) = B(T, - Tf = 5(1/ J, - I/Jf . (18) 

Using this simple form, without including any next-to-leading terms, we are able 
to fit all points of Table |l] with B = 0.9670(5) and the exponent P = 0.3785(6), in 
agreement with the high-precision zero-field determination in [|T^. 

As the critical point is reached the ii-dependence of the magnetization should 
change to satisfy critical scaling. We thus fit the data from the largest lattice sizes 
at Tc to the form 

M{T,,H) = d,H'/'. (19) 

As can be seen in Fig. ^a a very good straight-line fit to the largest-L results is 
possible. The smaller lattices show however definite finite-size effects. We find the 
exponent 6 = 4.86(1), in agreement with [0, and in addition we obtain the critical 
amplitude dc = 0.715(1). In Fig. ^ the magnetization at Tc is compared with the 
finite-size-scaling prediction 



(20) 



using the critical exponents of Ref. fl^. We observe no corrections to scaling, even 
at higher ii-values. The scaling function Qm is universal. In order to be consistent 
with Eq. ([191) for large z = HL^^^'^, i.e. for finite small H and large L, it must 
behave as 

QMiz) = d,z^". (21) 

This offers a second way to determine the critical amplitude dc-, this time exploiting 
also the data of the smaller lattices. From a fit in the z-range 20 — 1000 we find the 
value dc = 0.713(1), which agrees with our first determination. 

In Fig. ^ we show an example (at J = 0.98) of the different behaviours of xt and 
Xl at low temperatures. As a test we compare the result for dM/dH (line) from the 
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M-fits in Table |I| to the XL-data. Though there are large finite-size effects for small 
L, the results for the highest L-values agree nicely with the expected behaviour. 
A similar test can be done for xt by showing also the result for M/H (line), as 
obtained from the measurements of the magnetization in Fig. ^ Here as well we 
find agreement for large L. 
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Figure 4: The magnetization at Tc. In the left figure (a), M is plotted vs. H^/\ the 



line is the fit ([l9|) , the solid part shows the range used for the fit. The right plot (b) 
shows the finite-size-scaling function Qm from Eq. 




Figure 5: The susceptibilities below Tc at J = 0.98. On the left xl is plotted vs. 
^-1/2 rjj-^g right plot shows xt as a function of H^^. The lines are explained in 
the text. 
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The if-dependencies of the susceptibihties at the critical point are the same. 
Their amphtudes differ however by a factor 6. Here and in the following we use the 
values 6 = 4.86 and (3 = 0.38 . From Eqs. (§) and (P) and from the magnetization 
at the critical point, Eq. ([T9|), we derive 

XL = {dJ5)H'/'-' and xt = d,H^/'-^ . (22) 

In the left part of Fig. ^ we compare xl and xt at T^.. The lines in the figure are 
calculated from Eq. (|22|) and the fit results of Eq. ([T9|). We see again consistency with 
the highest-L data. The right part of Fig. ^ shows two examples of xl and xt for 
high temperatures. Both susceptibilities converge to one value x{T) for H 0, since 
no spontaneous symmetry breaking occurs for T > Tc. At the higher temperature 
corresponding to J = 0.80 the two susceptibilities are essentially constant (i.e. the 
magnetization is linear in H) and equal for a large range in H. However, at J = 0.90 
(that is closer to T^) and finite if- values xt is larger than xl- 
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Figure 6: The susceptibilities at Tc (left) vs. H^^^ ^ and above Tc (right) vs. H at 
J = 0.90 (squares) and J = 0.80 (circles). The lines are explained in the text. 



4 The Scaling Function 

The scaling function for the three-dimensional 0(4) model was determined from 
a fit of Monte Carlo data in [Q. Our goal here is to describe our data using the 
perturbative form of the equation of state as discussed in Section ^ but with non- 
perturbative coefficients, determined from a fit of the data. The original form of the 
function y = f{x) from [0] is not suitable for such a fit. We thus consider Eq. (|16|), 
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which is written as a simple series expansion in y. We do not expect this form to 
describe the data for all x and y, yet looking at Fig. [|, we hope to cover a significant 
portion of the phase diagram for small x. The idea is to interpolate this result with 



a fit to the large- a; form (17). 



Our fits are shown, together with our data, in Figs. |^ and ^. We have considered 
data from our largest lattices, for inverse temperatures 0.9 < J < 1.0 and magnetic 
fields H < 0.01. The normalization constants Hq and Tq, obtained from Eq. (JT^ ) 
and our fits in Section ^ are given by 

Ho = 5.08(3) , To = 1.093(2) . (23) 

We have performed a fit using the three leading terms in (|16D for small y 

xi{y) + 1 = {Zi + li)y + Z2y^'^ + rfs y'^' • (24) 

This form was fitted in the interval —1 < x ^ 1.5, giving 

0^ + 4 = 0.345(12) , Z2 = 0.6744(73) , d2 = -0.0232(49) . (25) 




Figure 7: The function y = f{x) from fits to data at small x (red line) and at large 
X (blue line). 



The fit describes all the data at T < T^, and also higher, up to a; ~ 5. This confirms 
that the expression (|I6|) is valid also away from x ~ — 1, as observed in Section |^. 
We note the small value of ^2- An attempt to include the next power of y leads to 
a coefficient that is zero within errors. We also see that our data are not sensitive 



to possible logarithmic corrections to Eq. ([Tq) as proposed in |jT5[. Our coefficients 
can be compared to those calculated perturbatively for = 4 in Ref. 

ci + 4 = 0.528 , C2 = 0.530 . (26) 
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For large x we have done a 2-parameter fit of the behaviour (|13), in the corre- 
sponding form for x in terms of y 

X2{y) = ay'/^ + . (27) 

Considering data points with y > 50 (corresponding to x > 15) we obtain 

a = 1.084(6) , b = -0.874(25) . (28) 

Expression ( PTj) is seen to describe the data for x ^ 2. We mention that a fit of our 
data to the leading term in Griffiths's condition, using x > 50, yields 7 = 1.45(1), 
which is in agreement with . 
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Figure 8: The scaling function fa 
together with our fit from Eq. (p9|). 



M/h^/^ from Eq. (|g). The dat shown 



The small- and large-x curves cover the whole range of values of x remarkably 
well. In fact, the two curves are approximately superimposed in the interval 2 < 
X < 8. We can therefore interpolate smoothly, for example by taking 



x{y) = xi{y) 



z/o + r 



+ X2iy) 



y 



y'o + y 



(29) 



at ?/o = 10, which corresponds to x ~ 4. Expression ( pPf ) is equivalent to the equation 
of state (|T3|) and to the scaling function fc in (|1^). In Fig. | we show a plot of fc 
obtained parametrically from x{y) in (^). The two variables in the plot are simply 
related to x and y by 



G 



M/h 



1/5 



-i/s 



t/h 



1/135 



xy 



-1/135 



(30) 



11 



We see a remarkable agreement of our data points with the form suggested by 
perturbation theory. With respect to the scahng function in our function is 
shghtly higher for large negative tjh^l^^ ^ due to our more complete elimination of 
finite-size effects. 



5 Summary and Conclusions 

We have shown that the Golstone singularities are clearly observable at low temper- 
atures, and also close to T^. In fact, we are able to use the observed Goldstone-effect 
behaviour to extrapolate our data to — and obtain the zero-field critical expo- 



nent /3 in good agreement with ll^. We remark that the same does not happen at 
high temperatures: we are not able to get the exponent 7 from extrapolations using 
the constant behaviour of the longitudinal susceptibility (or the linear behaviour of 
the magnetization), since this behaviour is masked close to for the fields H we 
have taken into account. At the same if's the H^l"^ behaviour is clearly present for 
all the T < Tf. consider, showing that the Goldstone effect is dominating the 
critical one, except at T^. 

A strong manifestation of the Goldstone behaviour had been conjectured pertur- 
batively based on the size of the coefficient C2 in the e-expansion of the equation 
of state. We have fitted the perturbative form in Section ^, finding a coefficient that 
is even larger than the perturbative one. 

The resulting curve for the equation of state describes all the data beautifully, 
and can be plotted parametrically for the scaling function. 

As a by-product of our work we have determined the critical exponent 5 = 4.86(1) 
by a fit of the magnetization at Tc to the critical scaling behaviour as a function of 
H. In addition we checked the finite-size-scaling prediction for M. It is remarkable 
that we observed in both cases no corrections to scaling. 



A similar investigation for the 0(2) model is currently being done |I8 
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